ERG and c-MYC regulate a critical gene network in BCR::ABL1-driven B cell acute lymphoblastic leukemia

Philadelphia chromosome–positive B cell acute lymphoblastic leukemia (B-ALL), characterized by the BCR::ABL1 fusion gene, remains a poor prognosis cancer needing new therapeutic approaches. Transcriptomic profiling identified up-regulation of oncogenic transcription factors ERG and c-MYC in BCR::ABL1 B-ALL with ERG and c-MYC required for BCR::ABL1 B-ALL in murine and human models. Profiling of ERG- and c-MYC–dependent gene expression and analysis of ChIP-seq data established ERG and c-MYC coordinate a regulatory network in BCR::ABL1 B-ALL that controls expression of genes involved in several biological processes. Prominent was control of ribosome biogenesis, including expression of RNA polymerase I (POL I) subunits, the importance of which was validated by inhibition of BCR::ABL1 cells by POL I inhibitors, including CX-5461, that prevents promoter recruitment and transcription initiation by POL I. Our results reveal an essential ERG- and c-MYC–dependent transcriptional network involved in regulation of metabolic and ribosome biogenesis pathways in BCR::ABL1 B-ALL, from which previously unidentified vulnerabilities and therapeutic targets may emerge.


INTRODUCTION
Acute lymphoblastic leukemia (ALL) is a highly aggressive cancer.Numerous genetic subtypes of ALL have been identified, including Philadelphia chromosome-positive B cell ALL (Ph + B-ALL) that comprises ~25 to 30% of all adult ALL.Ph + B-ALL is characterized by the t (9;22) chromosomal translocation that generates the BCR::ABL1 fusion gene resulting in abnormal tyrosine kinase signaling in B lymphoid progenitors that drives leukemia development.Ph + B-ALL remains a poor-prognosis cancer (1-3) despite the improved outcomes achieved by targeting the driver BCR::ABL1 fusion protein with tyrosine kinase inhibitors (TKIs) (4).More than 50 to 60% of patients relapse (5)(6)(7), including those receiving newer immunotherapeutic approaches.A better understanding of disease pathogenesis would assist rational development of new targeted therapeutic approaches.
While additional genomic lesions in BCR::ABL1 B-ALL have been identified, these are largely loss-of-function alleles, including the B cell transcription factor genes IKZF1 (IKAROS), PAX5, and EBF1, that are clinically associated with poorer disease outcomes (8)(9)(10)(11).Loss of these factors in preclinical models has been associated with accelerated leukemia expansion that is proposed to occur via loss of their metabolic gatekeeper function (12), suggesting that these genes and members of their gene regulatory networks are unlikely to represent viable molecular candidates as direct therapeutic targeting of loss-of-function alterations is problematic.
We sought to specifically identify gain of function and thus potentially targetable alterations that contribute to Ph + B-ALL pathogenesis, including those that may not be identified through diagnostic genomic assays.To do this, we used a murine model of BCR::ABL1 B-ALL to identify transcription factors and gene regulatory networks required for leukemia initiation and progression.Transcriptional profiling revealed up-regulation of ERG and c-MYC in murine and human BCR::ABL1 B-ALL, and, in murine and human models, ERG and c-MYC were required for BCR::ABL1 B-ALL.Profiling of ERG-and c-MYC-dependent gene expression and analysis of chromatin immunoprecipitation sequencing (ChIP-seq) data established that ERG and c-MYC coordinate a regulatory network in BCR::ABL1 B-ALL directly controlling expression of genes enriched for involvement in several cellular processes.Prominent was control of ribosome biogenesis, including direct regulation of RNA polymerase I (POL I) subunits, the role of which was validated by inhibition of BCR::ABL1 cells and other genomic subtypes by POL I inhibitors.

ERG, c-MYC, and downstream gene targets are highly expressed in Ph + B-ALL
We initially determined the transcriptional changes associated with leukemic transformation in BCR::ABL1 B-ALL in the P190 transgenic mouse model that expresses a human BCR::ABL1 transgene under the control of the metallothionein promoter (13).This transgenic mouse develops a highly penetrant acute pre-B cell (precursor B cell) lymphoblastic leukemia that closely phenocopies human disease.The murine model is characterized by expansion of pre-B cells in the bone marrow, typically between 5 and 8 weeks of age, and a B-lineage differentiation block resulting in a deficiency of maturing immunoglobulin M (IgM) + IgD + B cells (Fig. 1A).
Comparison of transcriptional changes in murine P190 leukemic pre-B cells to publicly available gene expression data from human BCR::ABL1 B-ALL (GSE5314) (22) revealed significant correlation with transcriptional changes in human B-ALL (Fig. 1, B and C).
Examination of transcription factors from the Kyoto Encyclopedia of Genes and Genomes deregulated in cancer gene set (ko05202) identified up-regulation of several oncogenic transcription factors involved in B cell differentiation, including the E26 Transformation-Specific (ETS) family transcription factor ERG (23), the β helixloop-helix transcription factor c-MYC (24) and its cooperative transcriptional partner MAX (25) in both murine and human BCR::ABL1 B-ALL (Fig. 1C).Increased protein levels of ERG and MYC were observed in P190 B-ALL cell lines compared to control pre-B cells (Fig. 1D).
Consistent with regulation by ERG and c-MYC of gene networks required for BCR::ABL1 B-ALL leukemogenesis, significant enrichment for the ETS, c-MYC, and MAX binding motifs was found in genes up-regulated in murine and human BCR::ABL1 B-ALL (Fig. 1E).Gene Ontology (GO) analysis of up-regulated genes in murine and human BCR::ABL1 B-ALL showed that the most significantly enriched gene sets were associated with ribosome biogenesis, as well as processes such as amide biosynthesis and nucleotide, tRNA, and amino acid metabolism (Fig. 1F and table S1).Up-regulation of ribosome biogenesis in BCR::ABL1 leukemia was confirmed by significantly increased transcription of the POL I-dependent 47S/45S precursor ribosomal RNA (rRNA) internal transcribed spacer in P190 leukemia cells compared to wild-type pre-B cells (fig.S1B).

ERG and c-MYC are critical for Ph + B-ALL development
To explore the functional significance of ERG and c-MYC in Ph + B-ALL, we first examined whether co-occurrence of BCR::ABL1 and ABL1 class fusions occurred with genomic variants of ERG and c-MYC in a cohort of human ALL (St.Jude, PeCan, accessed 21 January 2022) (26).While loss-of-function variants and copy number loss of ERG and, to a smaller extent, c-MYC were identified in this cohort of B-ALL, these did not co-occur with BCR::ABL1 or ABL1 class fusions (fig.S1C S1).
that, while gene dependency for ERG and c-MYC could be variably identified for several genomic subtypes in B-ALL cell lines, no BCR::ABL1 class cell lines had been assessed (fig.S1D).
To directly address the requirement for ERG and c-MYC in BCR::ABL1 B-ALL, we generated P190 mice in which the Erg or c-Myc genes were deleted from the Common lymphoid progenitor (CLP) stage of lymphopoiesis using a Rag1-Cre conditional knockout approach (28).Deletion of a single allele of Erg (P190 T/+ ;Rag1 Cre T/+ ;Erg ∆/+ ) was sufficient to prevent the development of P190 B-ALL in this model.Deletion of one c-Myc allele (P190 T/+ ;Rag1Cre T/+ ; c-Myc Δ/+ ) also significantly delayed leukemia development, and this delay was more pronounced in the absence of both c-Myc alleles (P190 T/+ ;Rag1Cre T/+ ;c-Myc Δ/Δ ) (Fig. 2A).At 5 weeks of age, pre-B cell numbers in P190 bone marrow were similar to wild-type mice, whereas, at 8 weeks of age, before overt symptoms of P190 disease, a significant proportion of P190 mice had developed an abnormal accumulation of pre-B cells (fig.S1E).In contrast, in P190 mice lacking one copy of Erg or c-Myc, no expansion of pre-B cells was seen, with pre-B cell numbers in 8-week-old P190 T/+ ;Rag1Cre T/+ ;Erg Δ/+ and P190 T/+ ;Rag1Cre T/+ ;c-Myc Δ/+ mice comparable to those seen in C57BL/6 controls (Fig. 2B).To determine the impact of loss of Erg or c-Myc alleles on clonal expansion during BCR::ABL1 leukemogenesis, immunoglobulin heavy-chain (Igh) gene clonotyping analysis was performed on bulk RNA-seq data obtained from primary pre-B cells.In P190 mice, leukaemogenesis was associated with dominant clones arising at 5 and 8 weeks of age, with one clone often becoming dominant in mice developing overt leukemia (Fig. 2C).In contrast, quantitative analysis of the 10 most frequent Igh clones revealed no significant dominant clonal expansion in P190 T/+ ;Rag1Cre T/+ ;Erg Δ/+ mice at 8 weeks of age.Similarly, dominant clonal expansion was not observed in P190 T/+ ;Rag1Cre T/+ ;c-Myc Δ/+ mice.Together, these data demonstrate high expression of ERG and c-MYC in BCR::ABL1 B-ALL compared to non-transformed pre-B cells and that these transcription factors are necessary for pre-B cell clonal expansion and subsequent leukemia development.

ERG and c-MYC contribute to BCR::ABL1 B-ALL maintenance
We next assessed the role of ERG and c-MYC in sustaining established BCR::ABL1 leukemia.We derived multiple independent cell lines from leukemic P190 mice carrying either floxed Erg (Erg fl/fl ) or c-Myc (c-Myc fl/fl ) alleles in addition to the CreERT2 transgene (29) to allow 4-hydroxy-tamoxifen (4-OHT)-dependent deletion of the floxed alleles (Fig. 3A, top).Upon 4-OHT treatment, reduced expression of ERG or c-MYC was observed in the respective cell lines, in which Erg or c-Myc was conditionally deleted (Fig. 3B).Dose-responsive inhibition of these ERG and c-MYC-deficient leukemic cell lines was observed, an effect not seen in in P190 T/+ ;CreER T/+ control leukemia cells (Fig. 3C).
We next investigated the requirement for ERG and c-MYC in BCR::ABL1 leukemia maintenance in vivo.Individual P190 T/+ ; CreER T/+ ;Erg fl/fl and P190 T/+ ;CreER T/+ ;c-Myc fl/fl leukemia cell lines were pretreated with either 4-OHT (pre4-OHT) resulting in significant loss of Erg or c-Myc expression (Fig. 3D) or vehicle control and then transplanted into irradiation-conditioned recipients.Mice that received vehicle-treated cells were then either given tamoxifen (TAM) at day 8 or left untreated (control group) (Fig. 3A, bottom).In recipients of cells that are ERG-or MYC-deficient, by either 4-OHT pretreatment or in vivo TAM administration, the time to ethical endpoint due to leukemia was prolonged relative to control mice (Fig. 3E).Splenomegaly was observed to be significantly reduced in cohorts that received c-MYC-deficient cells (Fig. 3F), while the proportion of donor cells in bone marrow was observed to be consistently high in all groups (Fig. 3G).Notably, leukemia that developed in mice transplanted with P190 T/+ ;CreER T/+ ;Erg fl/fl and P190 T/+ ;CreER T/+ ;c-Myc fl/fl cell lines following 4-OHT treatment or TAM treatment in vivo demonstrated incomplete and variable reduction in ERG or c-MYC expression in the diseased bone marrow cells (fig.S2A), suggesting in vivo expansion of leukemia cells that had escaped efficient Cre-mediated gene recombination had occurred.
To confirm and extend our observations to human BCR::ABL1 B-ALL cells, guide RNAs directed against human ERG or c-MYC were expressed via a doxycycline (Dox)-inducible lentiviral vector co-expressing green fluorescent protein, resulting in reduced ERG or c-MYC protein levels in respective Cas9-expressing human BV173 cells (Fig. 4A).Reduction of ERG or c-MYC expression by two independent guide RNAs for each gene resulted in a distinct competitive proliferative disadvantage in BV173 cells compared to empty vector controls (Fig. 4B).Similar results were obtained in a second BCR::ABL1 B-ALL human cell line, SupB15 (fig.S2, B and  C).Last, reduction of ERG or c-MYC either before transplantation (preDox) or treatment with Dox in vivo resulted in delayed BV173 tumor growth in transplanted mice compared to untreated controls or mice transplanted with empty vector expressing BV173 cells (control) (Fig. 4C).
In addition to BCR::ABL1 B-ALL, prominent ERG and c-MYC expression is also found in other B-ALL subtypes (fig.S3A) and pre-B-ALL human cell lines, including Nalm6 and RS4:11 (fig.S3B), where binding of ERG and c-MYC to the POLR1B promoter can be observed (fig.S3C).As observed in BCR::ABL1 B-ALL cell lines, genetic reduction of ERG or c-MYC expression in Nalm6 and RS4:11 cells resulted in a competitive proliferative disadvantage compared to empty vector controls (fig.S3, D to G).

Identification of a gene network regulated by ERG and c-MYC in BCR::ABL1 B-ALL
To define the gene networks regulated by ERG and c-MYC that facilitate leukaemogenesis, we first examined the gene expression changes upon deletion of either Erg or c-Myc in established P190 T/+ ; CreERT2 T/+ ;Erg fl/fl and P190 T/+ ;CreERT2 T/+ ;c-Myc fl/fl cell lines (table S2).Notably, there was overlap of differentially expressed genes upon Erg or c-Myc deletion in genetically independent cell lines.Of particular interest were genes down-regulated with both Erg or c-Myc deletion (Fig. 5A), as these genes may form part of a transcriptional gene network regulated by ERG and c-MYC, mediating the functional roles for these transcription factors in BCR::ABL1 B-ALL during leukaemogenesis and in leukemia maintenance.To define genes directly bound and transcriptionally regulated by ERG and c-MYC, publicly available pro-B cell ChIP-seq datasets for ERG (GSM3895108) and c-MYC (GSM1234475) were examined.An initial genome-wide motif analysis of the ERG ChIP-seq dataset identified enriched representation of the ERG motif at ERG binding sites as expected (Fig. 5B).In addition, enrichment of the β helixloop-helix binding motif recognized by c-MYC (30) and the YY1 cohesin motif were also identified within 50 base pairs (bp) of ERG binding sites.For the complimentary analysis for c-MYC binding, the most highly enriched motifs within 50 bp of c-MYC-bound loci were not only motifs for the c-MYC binding partner MAX (31) as well as c-MYC, as expected, but also the ERG motif.We next focused on analysis of genome wide ERG and c-MYC binding to regions within 5 kb of a transcriptional start site (TSS).This analysis demonstrated overlap of ERG and c-MYC binding to a subset of genomic loci (Fig. 5C) and suggested that c-MYC and ERG co-localization at specific genomic loci may result in co-regulation of specific genes as part of a transcriptional network.Notably, no direct interaction between ERG and c-MYC was identified by co-immunoprecipitation (fig.S4A).
To explore this hypothesis further in BCR::ABL1 B-ALL, differentially expressed genes that were found to be down-regulated upon both Erg and c-Myc deletion (Fig. 5A) were compared to genomic loci bound by ERG and c-MYC to define target genes of the transcriptional network directly co-regulated by c-MYC and ERG (table S3).GO analysis was then undertaken to define the biological functions of these genes.This identified that the network was significantly enriched for genes involved in ribosome biogenesis as well as for genes involved in several metabolic processes (Fig. 5D).These findings were congruent with the GO analysis of up-regulated genes in human and murine BCR::ABL1 B-ALL (Fig. 1F), thereby establishing the key roles for c-MYC and ERG in regulating these leukemiaassociated changes.These data also identified direct regulation by ERG and c-MYC of genes involved in ribosome biogenesis such as the nucleolar proteins nucleophosmin (Npm1) (32) and fibrillarin (Fbl) (33) and subunits of the POL I complex (Polr1b and Polr1c) responsible for transcribing rRNA genes within the nucleoli (34) (Fig. 5E and table S3).This was in contrast to examples of genes specifically regulated by ERG or c-MYC (fig.S4B).
As molecular validation of this analysis, we targeted POL I to disrupt ribosome biogenesis in BCR::ABL1 B-ALL using the firstin-class small-molecule inhibitor of POL I transcription, CX-5461 (Pidnarulex).POL I transcription occurs in the nucleoli and produces the precursor rRNA (pre-rRNA) containing the sequences of the 18S, 5.8S, and 28S mature rRNA components of the ribosome (35,36).We confirmed on-target CX-5461-mediated inhibition of POL I transcription by measuring 47S/45S pre-rRNA levels in BCR::ABL1 B-ALL cell lines (Fig. 6A).As inhibition of POL I transcription results in nucleolar disruption and the induction of nucleolar stress response, we assessed the effect of CX-5461 on nucleolar integrity using immunofluorescence staining for the nucleolar protein fibrillarin (FBL), a small nucleolar ribonucleoprotein that directs the methylation and processing of pre-rRNAs (37)(38)(39).In murine and human BCR::ABL1 B-ALL cell lines, CX-5461 treatment resulted in reduced FBL staining  S3) with Polr1b, Polr1c, Fbl, and Npm1 nucleolar genes highlighted.(E) RnA-seq and chiP-seq tracks of Polr1b, Fbl and Npm1 gene loci showing transcriptional changes associated with erg (P190 T/+ ;CreER T/+ ;Erg Δ/Δ ) and c-Myc (P190 T/+ ;CreER T/+ ;Myc Δ/Δ ) deletion compared to P190 T/+ ;CreER T/+ ;Erg fl/fl and P190 T/+ ;CreER T/+ ;Myc fl/fl controls, and erg and c-Myc binding to promoter regions (highlighted in blue) defined by h3K4me3 (GSM2255547) and h3K27Ac marks (GSM2255552), with tracks for independent chiP-seq undertaken in the P190 BCR::ABL1 cell line for eRG, h3K4me3, and h3K27Ac shown.
that was localized in smaller punctate nucleolar domains as well as diffusely in the nucleoplasm, indicating significant nucleolar disruption (fig.S5, A and B).Consistent with the established action of CX-5461 in activating the nucleolar stress response (39,40), CX-5461 treatment was associated with up-regulation of the tumor suppressor TP53 (fig.S5C).CX-5461 treatment of murine and human BCR::ABL1 leukemia cells caused dose-dependent growth inhibition at nanomolar drug concentrations (Fig. 6B), which has previously been shown to be due to induction of apoptosis in TP53 replete leukemia cells (37).Sensitivity to inhibition of POL I transcription was confirmed with two other POL I inhibitors, actinomycin D (41) and BMH-21 (42) (Fig. 6C).CX-5461 also demonstrated activity across other B-ALL genomic subtypes with dose-dependent growth inhibition of the DUX4::IGH (Nalm6) and KMT2A::AF4 (RS4:11) B-ALL cell lines (fig.S5D), an observation that could also be extended with other POL I inhibitors (Fig. 6C).

DISCUSSION
Ph + B-ALL remains a poor-prognosis genetic subtype of leukemia.The defining t(9;22) translocation generates the oncogenic BCR::ABL1 fusion protein that is the molecular driver of this leukemia and the primary target for small-molecule kinase inhibitors.However, despite the advent of TKI therapy combined with conventional intensive chemotherapy (4), leukemia relapse often occurs from preexisting clones carrying mutations in the tyrosine-kinase domain of the BCR::ABL1 fusion gene (43) that render targeted kinase inhibition ineffective.There is a need to develop more effective therapeutic strategies to treat this disease, including identification of molecular targets other than BCR::ABL1 that can bypass resistance to TKI inhibitors and other potential mechanisms of chemotherapy resistance.
Transcriptional profiling of human and murine BCR::ABL1 B-ALL identified members of two sequence-specific transcription factor families whose expression was high in Ph + B-ALL: ERG, a member of the ETS family, and the β helix-loop-helix family member c-MYC and its heterologous binding partner MAX (44).Notably, both ERG and c-MYC have also been shown to have critical roles in B cell development (23,24).murine and human BCR::ABL1 B-ALL cell lines, reduction of ERG or c-MYC levels constrained leukemia cell expansion and delayed leukemia development in transplanted mice.Examination of publicly available gene expression data suggests that ERG and c-MYC are co-expressed in other B-ALL subtypes in addition to BCR::ABL1 B-ALL.Our observation that genetic reduction of ERG or c-MYC resulted in proliferative disadvantage in two other pre-B-ALL genomic subtypes provides preliminary evidence that our observations may extend more broadly to other forms of B-ALL, an observation that was also noted with variable gene dependency for ERG and c-MYC in B-ALL cell lines of different genomic subtypes in the Cancer Dependency Map initiative (27).
Genome-wide analysis of ERG-and c-MYC-dependent transcriptional changes in human and mouse BCR::ABL1 B-ALL cell lines, combined with analysis of ChIP-seq datasets, defined a cooperative ERG-and c-MYC-dependent transcriptional network that included direct binding and transcriptional regulation of genes involved in metabolic pathways and ribosome biogenesis.Enrichment of genes involved in ribosome biogenesis, including components of the POL 1 transcription complex, allowed a molecular validation of the network via inhibition of POL 1 transcription using CX-5461, actinomycin D, and BHM-21, each of which caused dose-dependent growth inhibition of BCR::ABL1 B-ALL cell lines at nanomolar drug concentrations.CX-5461 has previously been shown to selectively induce apoptosis in c-MYC driven leukemia cells (37).While an oncogenic role for c-MYC has been proposed via transcriptional amplification (45,46), prior evidence for c-MYC as a master regulator of ribosome biogenesis (47,48) and the observation that c-MYC-driven cancers are associated with hyperactivated POL I transcription (35,36) are particularly supported by our data, which now extends a similar role for c-MYC to BCR::ABL1 B-ALL.While, in hematological disease, ERG has established oncogenic roles in myeloid malignancy (49)(50)(51)(52)(53) and T cell ALL (54,55), our finding of its cooperative role with c-MYC was unexpected and provides an important insight into the contribution of these two transcription factors in the biology of BCR::ABL1 B-ALL.The transcriptional regulation by ERG and c-MYC of genes involved with metabolic pathways and ribosome biogenesis (56) may have broader implications for other malignancies, in which ERG is deregulated.
In summary, we have established essential roles for c-MYC and ERG in BCR::ABL1-driven B-ALL and defined an ERG-and c-MYC-dependent transcriptional network involved in regulation of metabolic processes and ribosome biogenesis in this disease.Together, our results validate an approach for defining essential transcriptional regulatory networks to elucidate important biological pathways in oncogenesis, among which previously unidentified vulnerabilities and therapeutic targets may emerge.

Mice
Unless otherwise described, mice were generated on a C57BL/6 background.Mice with a conditional Erg knockout allele (Erg fl ) were generated as previously described (23).Mice carrying a conditional c-Myc knockout allele (24) (Myc fl ) were obtained from the Jackson Laboratory (Myc tm2Fwa ).These mice were interbred with either Rag1Cre mice (28), in which Cre recombinase is expressed during lymphopoiesis from the CLP stage (57) or CreERT2 mice (29), in which the expression of the Cre recombinase can be initiated by TAM treatment to generate Rag1Cre T/+ ;Erg fl/fl , Rag1Cre T/+ ;Myc fl/fl , CreERT2 T/+ ;Erg fl/fl , CreERT2 T/+ ;Myc fl/fl , and wild-type control littermates.Subsequently these were crossed to P190 transgenic mice (13).Mice were co-housed in a barrier facility and analyzed from 6 to 18 weeks of age.Male and female mice were used.The primers used for genotyping are provided in table S4.NOD.Cg-Prkdc scid IL2rg tm1Wjl / Szj (NSG) female mice were obtained from the Jackson Laboratory and co-housed in individually ventilated cages in a specific pathogenfree facility.Experimental procedures were approved by The Walter and Eliza Hall Institute of Medical Research Animal Ethics Committee.

In vitro cell assays
To measure adenosine 5′-triphosphate levels after 4-OHT, CX5461, actinomycin D, or BMH-21 treatment, 2500 cells were seeded per well in triplicates in the indicated drug concentrations and analyzed after 48 to 96 hours of culture as indicated by the CellTiter-Glo Luminescence Assay (Promega).For in vitro competition cultures, cells transduced with a Dox-inducible fluorescently labeled short hairpin RNA/guide RNA vector were co-cultured with uninfected parental cells, split into two aliquots and treated with Dox (100 ng/ ml) or dimethyl sulfoxide (DMSO).The relative frequency of transduced and non-transduced cells was measured over time by flow cytometry.

Flow cytometry
For flow cytometric analysis, single-cell suspensions were prepared in balanced salt solution (0.15 M NaCl, 4 mM KCl, 2 mM CaCl 2 , 1 mM MgSO 4 , 1 mM KH 2 PO 4 , 0.8 mM K 2 HPO 4 , and 15 nM Hepes supplemented with 2% bovine calf serum).Cells were washed, stained with fluorophore-conjugated antibodies (see table S5), and analyzed or sorted on a BD LSRFortessa or BD FACSAria III, respectively.Dead cells were excluded by staining with FluoroGold (AAT Bioquest), and data analysis was performed using FlowJo 10.4 (Becton Dickinson).To determine total cell numbers, an aliquot of the single-cell suspension was mixed with a defined number of allophycocyanin (APC + ) beads.The ratio of cells/bead was used to determine the total cell count per femur.

Immunoblots and co-immunoprecipitation
For Western blotting, cells were lysed in KALB lysis buffer [1% Triton X-100, 150 mM NaCl, 50 mM tris-HCl (pH 7.4), 1 mM EDTA, 1 mM phenylmethylsulfonyl fluoride, 2 mM Na 3 VO 4 , 10 mM NaF, and complete protease inhibitors (Roche)].Proteins were separated in 4 to 12% Bis-Tris NuPAGE protein gels (Invitrogen) under reducing conditions, transferred onto a Immobilon-P membrane (Millipore), and immunoblotted with primary antibodies (see table S5), followed by incubation with secondary horseradish peroxidase-conjugated antibodies and visualization by enhanced chemiluminescence.For co-immunoprecipitations, cells were lysed as described above.Cell lysates were incubated with 0.25 μg of specific antibody or IgG control for 2 hours, Protein A slurry was added, and the mix was incubated for >3 hours before Protein A beads were washed three times in KALB lysis buffer and protein was eluted and then analyzed by Western blot.

Real-time PCR analysis
rRNA transcription by POL I in mouse and human BCR::ABL1 leukemia cells was assessed at the 47S/45S rRNA ITS and 47S/45S rRNA 5′ETS genes relative to β2m and VIM genes, respectively (see table S4).RNA was reverse-transcribed to cDNA using random hexamer primers (Promega, Madison, USA) and Superscript III reverse transcriptase (Invitrogen) according to the manufacturer's instructions.Quantitative polymerase chain reaction (PCR) was performed using SYBR green reagents (Applied Biosystems, USA) on the ViiA 7 real-time PCR system (Thermo Fisher Scientific, USA), and relative expression of rRNA genes after treatment with CX-5461 compared to housekeeping genes was expressed relative to DMSO-treated cells using the 2 −ΔΔCT method (59).

RNA sequencing
Total RNA was extracted using the RNeasy Plus minikit (Qiagen) from murine BCR::ABL1 B-ALL cell lines and B220 + selected pre-B cells.RNA (500 ng) was used to generate cDNA libraries using TrueSeq Stranded mRNA kits (Illumina).Sequencing was performed on a Hi-Seq2500 or NovaSeq sequencing system (Illumina) to generate 100 bp single-end reads.Pseudo-biological replicates were sequenced for each BCR::ABL1 B-ALL cell line.Reads were aligned to the GRCm38/mm10 build of the Mus musculus genome using Rsubread align function, and read counts were summarized at the gene level (60).Genes were filtered as non-expressed if they were assigned 0.5 counts per million mapped reads (CPM) in fewer than two libraries.Counts were transformed to log 2 -CPM and the meanvariance relationship estimated using the voom function in limma (61).Library sizes were trimmed mean on M-values (TMM)normalized, and differential expression was assessed using quasilikelihood F tests (62).Genes were called differentially expressed if they achieved a false discovery rate of 0.05.For plotting purposes, counts were converted to reads per kilobase of transcript per million mapped reads (RPKM) using the rpkm function in limma.These data have been deposited in Gene Expression Omnibus database (accession number GSE213791).For clonotype analysis, fastq files from pre-B cells from individual mice or pooled from n = 3 mice for samples taken at 5 weeks and primary leukemia were analyzed using the MiXCR software package (3.0.6) (63).The frequency of the 10 most prevalent clonotypes was normalized to frequency per mouse.

Gene network analysis
All ERG and c-MYC ChIP-seq peaks mapping to differentially expressed genes in both P190 T/+ ;CreERT2 T/+ ;Erg Δ/Δ and P190 T/+ ;CreERT2 T/+ ; c-Myc Δ/Δ cell lines within 10 kb of the TSS were identified.Peaks inside the gene body were annotated as "proximal targets, " peaks overlapping the TSS were labeled as promoter regulated targets, peaks less than 3 kb upstream or downstream of the TSS were labeled as putative promoter regulated targets, and peaks more than 3 kb upstream or downstream TSS were labeled as putative distal targets (see table S3).GO annotation of differentially expressed genes was performed and underwent expert manual curation.The network was constructed using igraph CRAN package (66) and exported to Cytoscape (67) for customization using RCy3 (68) R/Bioconductor package.

Immunofluorescence analysis
Suspension cells were fixed in a 4% paraformaldehyde solution for 10 min and then cytospun onto pre-charged Super Frost Plus Slides (Menzel Gläser) using a double cytology funnel.Slides were permeabilized for 10 min in ice-cold phosphate-buffered saline (PBS; 0.1 M, pH 7.4) containing 0.3% Triton X-100.Following permeabilization, cells were washed three times with PBS (5 min each with gentle rocking on a lab shaker) and blocked in PBS consisting of 5% goat serum and 0.3% Triton X-100 for 30 min at room temperature.Next, slides were incubated with primary antibodies (see table S5) in 1% bovine serum albumin (BSA)/PBS for 1 hour at 37°C in a humidified chamber.Following staining, slides were washed three times with PBS (5 min each with gentle rocking) and then subsequently incubated with secondary antibodies (see table S5).Secondary antibodies were diluted 1:600 in 1% BSA/PBS, and slides were stained for 1 hour at 37°C.Last, cells were washed in PBS and counterstained for 10 min with 4,6-diamidino-2-phenylindole (1 μg/ml; Sigma-Aldrich) before being mounted with glass coverslips using VECTASHIELD Antifade Mounting Medium (Vector Laboratories).
Fluorescent confocal images were acquired using a Nikon C2 laser scanning confocal microscope system (Nikon, Melville, NY) equipped with a 60× oil immersion objective and NIS-Elements software (Nikon, Melville, NY) for acquisition of the images.A maximal intensity projection of a Z-stack was than generated using the software program ImageJ (1.47v, National Institutes of Health).Images were analyzed using CellProfiler version 3.1.9(Broad Institute) using the same manually set parameters and thresholds.For statistical analysis, mean signal intensity data values were normalized to the median of each respective vehicle control and the data plotted using GraphPad Prism software (version 7) performing a two-sided unpaired Mann-Whitney test where appropriate.

Supplementary Materials
This PDF file includes: Figs.S1 to S5 legends for tables S1 to S3 tables S4 and S5 legend for data S1 Other Supplementary Material for this manuscript includes the following: tables S1 to S3 data S1

Fig. 1 .
Fig. 1.Gene expression analysis in BCR::ABL1 B-ALL.(A) Representative flow cytometry plots of bone marrow B-lymphoid populations from c57Bl/6 wild-type (Wt) control mice and leukemic P190 mice with the percentage of cells in each gated population indicated.the immunoglobulin d (igd)/igM profile is from B220 + gated cells; the cd24/cd43 profile from B220 + igd − igM − gated cells.(B) Barcode enrichment plots of up-regulated (red) and down-regulated (blue) genes in human Ph + B-All (GSe5314, n = 37) compared to murine P190 B-All gene expression changes (n = 5).P values were calculated by ROASt gene set testing.(C) Scatter plot of log 2 fold change (logFc) in expression of orthologous genes in human Ph + B-All (GSe5314) and murine P190 B-All compared to normal pre-B cells.transcription factors from Kyoto encyclopedia of Genes and Genomes transcriptional deregulation in cancer (ko05202) and the quadrant of up-regulated genes in both human and murine BCR::ABL1 B-All are highlighted in red.(D) Western blot analysis of Abl, erg, and c-Myc in pre-B cells, compared to three independent P190 T/+ cell lines.Actin serves as loading control.(E) in silico motif discovery of up-regulated genes in human and murine BCR::ABL1 B-All [enrichr, https://maayanlab.cloud/enrichr/; with UcSc transcription factor encOde motif discovery (72) and encOde consensus target Genes and cheA transcription Factor chiP-X (73)].(F) Gene Ontology (GO) analysis of up-regulated genes in human and murine BCR::ABL1 B-All (Metascape, https://metascape.org;see also tableS1).